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Abstract 

We study the time evolution of configurations in the form of two parallel domain walls 
moving towards each other in a supersymmetric field model. The configurations involved are not 
BPS-saturated. It is found that for such collisions there exists some critical value v cr ~ 0.9120 
of the initial velocity t>j of the walls. At Vi < v cr we observed reflection, that was not followed 
by change of vacuum states sequence. In collisions with Vi > v cr the sequence of vacuum 
states changes. The results of the numerical simulations are in agreement with "potential" 
consideration. 
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The dynamic properties of domain walls in supersymmetric theories have attracted some 
attention recently 0J - fj. Depending on particular form of the superpotential being chosen, 
one obtains different sets of supersymmetric vacua and different structure of domain wall con- 
figurations interpolating between them. We restrict ourselves by consideration of the theory 
described by the superpotential 

W($,X) = — $ - -A$ 3 -a$X 2 , (1) 
A 3 

where m is a mass parameter and a and A are coupling constants. We assume that a and A are 
real and positive. The Lagrangian for the real parts of the scalar fields is given for this theory 
by the expression 

L = {d<pf + {d X f A0 2 - «x 2 ) 2 - 4a W ■ (2) 

The potential term of Eq. (2) has four degenerate vacuum states, shown in Fig. 1. This 
theory possess a wide class of domain walls, which link different vacua. Some of them satisfy 
first order differential equations analogous to the Bogomol'nyi-Prasad-Sommerfeld (BPS) equa- 
tions 0. The dynamic properties of the BPS configurations for this model were intensively 
studied recently |IJ - 0. This our work is devoted to the so-called non-BPS domain walls, i.e. 
configurations which link different vacua of the theory, but do not satisfy BPS equations. 
It is convenient to work with dimensionless field variables / and h, defined as 

m . m . 

= -rf , X = —?^ h . 
A v Aa 



The Lagrangian (2) yields the following equations of motion for fields / and h: 

P 

hi - V 2 h - -h{\ -f- h') + j 2 hf = . (3) 

Here p = A/a, m = 1. It was shown (see, e.g., Ref. ||) that for the case p = 4 the field 
equations (3) possess an "elementary" walls, connecting vacua 3 and 2, 2 and 4. Their form 
may be obtained analytically ||: 



/32(«) = ^(l + taoh|), M2)=y^(l-tanh£) ; (4) 



fu(z) = - (l - tanh , h 2A (z) =-J-(l + tanh , (5) 

here z is a space coordinate orthogonal to the walls. It is easy to see, that the rest energy of 
these 3 — > 2 and 2^4 walls equals E Q = 4/3. 
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Consider a non-BPS ansatz configuration 3 — > 2 — > 4 constructed from two elementary 
domain walls 3 — ► 2 and 2 — > 4 located at z = —z and z = +2 respectively. Let us take their 
simple superposition in the form 

/324(Z, ^o) = /32(^ + 2b) + /24^ - Z ) - 1 , 

^324(2, ^0) = ^32(2 + z ) + h 24 (z - z ) . (6) 

Note, that from the system (3) the special "diagonal" solution 3 — > 4 can be easily found by 
substituting / = 0: 

f 34 (z) = 0, h 34 (z) = -tanh^ . (7) 

The energy of such configuration is E 34 = 16/3. 

To get ^o-dependence of the energy of configuration (6) we have to insert (6) into Hamilto- 
nian of the model. As a result we obtain 

£324(20) = 2E + AE 324 (z ) , (8) 

where 

-4-rvi _ 

AE - I dz ^ C ^ 32 C ^ 24 ' n-dh 32 dh 24 , M j, 2 i_2 \2 , 4 



/* , «/32«/24 „ dll 32 dll 24 , 2 7,2 n2 , 4 <-2 i2 

324 -y ^ 2 ^717 + 2P ^7^7 + (1_/324 _/l324) + /324^324 

— oo 

— (1 — ^32 — ^32) — "^32^32 — (1 — ^24 ~~ ^24) — "^24^24 



(9) 



Here / 32 = / 3 2(^ + z ), h 32 = h 32 (z + z ), f 24 = h±{ z ~ ^o), ^24 = ^24(2 - 20) ■ We calculated 
the z -dependence of AE 324 numerically, see Fig. 2 (solid curve). At the limit of large zq the 
configuration (6) looks like two isolated walls 3^2 and 2 — > 4. Therefore their total energy 
equals 2E , and AE 324 « 0. As it is seen from Fig. 2, energy AE 324 increases with decreasing 
zq. At z = AE 324 (0) m 3.119. It corresponds to £324(0) = 2E + AE 324 (0) m 5.786. Note, 
that E 324 (0) is larger than E 34 = 16/3 w 5.333. The energy of configuration 3 — > 2 — > 4 (6) has 
its absolute maximum at z m —0.37 when (AE 324 ) max rs 3.202. At large negative zq AE 324 (zq) 
has asymptotic value about 2.274. In the range z < configuration (6) actually has the shape 
of the 3 — > 1 — > 4 type, see Fig. 3. 

It is clear, that we can construct an ansatz configuration 3 — > 1 — > 4 in analogy to (6): 

/314(2, z ) = / 3 i(z + z Q ) + / 14 (z - z ) + 1 , 

^314(2, Zq) = h 31 (z + Zq) + h 14 (z - Zq) , (10) 

where 



/si(«) = ~(l + tanh|), M«) = ^(l-tanh|) ; (11) 
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/ u (2) = --(l-tanh|), h u {z) = -J-(l+taah^ . (12) 

The energy of 3 — > 1 and 1 — > 4 walls is exactly the same as of 3 — > 2 or 2 — ► 4. Hence, the 
z - dependence of the energy of configuration (10) will be 

£ 314 (2 ) = 2E o + &E 3U (z ) , (13) 

with "potential" AE^^Zq) analogous to AE 32 i(z ), Eq. (9). The shape of AE 3U (z ) is the 
same as of AE 3 2a(zq). As it was already mentioned, at zq < ansatz (6) has the shape 
of 3 — > 1 — > 4 type indeed. Obviously, configuration (10) at negative z has the shape of 
3 — > 2 — > 4 type, see Fig. 3. Notice, that h^^z, —z ) ^ h 3 u(z, z ). If we would like to compare 
energies of configurations (6) or (10) that belong to one of these two types, we have to place 
curves AE 32 4(z ) and Ai? 314 (— zq) (or curves AE 31i (z Q ) and AE 32A {— z )) in the same plot. 
Fig. 2 is constructed just in this way. 

We solved field equations (3) numerically with initial conditions in the form of (6), where 
3^2 and 2^4 walls located at some initial distance 2z ^> 1 and are moving towards each 
other with some initial velocity V{. Depending on the initial velocity we observed different types 
of evolution. If Vi is less than some critical value v™ m , walls 3 — > 2 and 2 — > 4 collide and then 
escape from each other to infinity. As a result we return to the configuration of the 3 — > 2 — > 4 
type. At initial velocities > v™ m the walls collide in a different way. The point is that after 
collision the configuration of the type 3 — ► 1 — > 4 appears. From these numerical simulations 
we found v™ m ~ 0.9120. The presence of different regimes in such collisions is a consequence 
of the fact that the energy of configuration 3 —>■ 2 — ► 4 is not degenerate with respect to the 
parameter z$. So, we have a kind of "potential" interaction between 3^2 and 2^4 domain 
walls. It is worth to mention here, that in the case of BPS-saturated (or simply BPS) walls |J 
there is no potential interaction. The latter property is a consequence of the degeneracy in 
energies of configurations with different interwall distances, analogous to our parameter Zq. 

Existence of the critical velocity can be understood in terms of the potential approach. 
From Fig. 2 it is seen, that if the initial kinetic energy of the walls 3 — > 2 and 2 — > 4 is smaller 
than AE* pa 3.119, then (inelastic) reflection may be expected. If the kinetic energy of the 
walls exceeds AE*, it is natural to expect that configurations of the type (6) with negative zq 
appear. But configuration (6) at negative z is of the type 3 — > 1 — >• 4 indeed, and from Fig. 2 
we see, that in this sector configurations (10) have smaller energy. Hence, configuration (6) at 
negative Zq transforms into (10). In further evolution the walls 3—^1 and 1 — >• 4 escape to 
infinity. 

It is also worth mentioning, that the initial configuration (6) with zq = and v j = looks 
like some excitation over the static solution (7). After emission of part of energy in the form 
of waves, the evolution of this initial configuration (zq = 0, V{ = 0) leads to formation of an 
excited kink of type (7) (wobbling kink). We were unable to get this wobbling kink solution 
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making numerical calculations of the equations of motion (3) with initial conditions (6) when 
either z or Vi was not equal to zero. 
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Figure captions 



Fig. 1. Locations of the vacuum states of the model. 

Fig. 2. The profile of the potential A£ , 32 4 versus zq (solid curve) and the profile of the 
potential A.E314 versus (—Zq) (dashed curve). 

Fig. 3. Profiles of f(z) (solid lines) and h(z) (dashed lines) for configurations 3 — ■> 2 — > 4 
and 3 -> 1 -»• 4 at z = ±10.0. 
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